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We propose a novel Hamiltonian engineering method in spin-based quantum information pro- 
cessing architectures that requires no local control but only relies on collective spin rotations and 
field gradients. The technique relies on coupling engineering via the time-domain construction of a 
Bragg grating and a weighting function that modulate the coupling strengths. As an example, we 
demonstrate how to engineering the optimal Hamiltonian for perfect quantum information trans- 
port between two separated nodes of a large spin network. We engineer a spin chain with optimal 
couplings from a large spin network, such as naturally occurring in crystals, while decoupling all 
unwanted interactions. We show that for realistic experimental parameters this Hamiltonian engi- 
neering method could be used for perfect quantum information transport at room-temperature. 
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Introduction - Controlling the evolution of complex 
quantum systems has emerged as an important area of 
research for its promising practical applications. Such 
control can be often reduced to the goal of Hamiltonian 
engineering [1] (also extended to reservoir engineering [2- 
4]). Hamiltonian engineering has been used to achieve a 
variety of tasks, including quantum computation [5], im- 
proved quantum metrology [6] and protection against de- 
coherence by dynamical decoupling filtering [7, 8]. The 
most important application is quantum simulation [9], 
as first proposed by Feynman [10]. The ultimate goal is 
to achieve a programmable universal quantum simulator, 
able to mimic the dynamics of any system. A possible 
strategy is to use a quantum computer and decompose 
the desired evolution into unitary gates [11, 12]. Alterna- 
tively, one can use Hamiltonian engineering by a Suzuki- 
Trotter factorization of the desired interaction into ex- 
perimentally achievable Hamiltonians [13, 14]. However, 
experimental implementations of these simulation meth- 
ods often require local quantum control, which is difficult 
to achieve in large quantum systems. 

Instead, here we present a general scheme for quan- 
tum simulation that employs only collective rotations 
of the qubits and field gradients - technology which is 
in everyday use e.g. in magnetic resonance imaging or 
in optical lattices [15]. We consider a qubit network 
with an internal Hamiltonian H lnt , for example due to 
the dipolar coupling between spins in a crystal lattice. 
The target Hamiltonian is then engineered from H m t 
in two steps: First, unwanted couplings are removed 
{decoupled) by evolution under the gradient field, in a 
manner equivalent to the construction of a time-domain 
Bragg grating Gij [7]. The grating acts like a sharp fil- 
ter, letting through only specific couplings [16]. Then, 
a time-domain weighting function F modulates the cou- 
pling strengths to match the target Hamiltonian. 

For concreteness, we show how to apply this filtered 
engineering method to the specific problem of generating 




FIG. 1. A complex spin network in a trigonal planar lattice 
consisting of a linear chain (red) and surrounding network 
(orange). The edges denote couplings between the spins. A 
linear magnetic field gradient is aligned in the direction of 
the chain, such that the gradient strength at spin j is Uj = 
jcj. The filter eliminates off-chain couplings as well as NNN 
couplings. 



an optimal Hamiltonian for quantum information trans- 
fer (QIT) between two separated nodes of a spin network. 
Linear arrays of spins have been proposed as quantum 
wires for this task [17] and engineering the coupling be- 
tween the spins can achieve perfect QIT [18]. We will 
finally analyze experimental requirements to implement 
the method in existing physical architectures. 

Hamiltonian engineering - The goals of filtered 
Hamiltonian engineering can be summarized as (i) the 
cancellation of unwanted couplings - often next-nearest 
neighbor (NNN) and other long range couplings - and 
(ii) the engineering of the remaining couplings - for ex- 
ample, nearest neighbor (NN) couplings - to match the 
desired coupling strengths. We achieve these goals by 
dynamically generating tunable and independent grating 
(Gij) and weighting (Fij) functions by means of collec- 
tive rotations and gradients. The first step is to create 
from Hi nt the Hamiltonian operator one wishes to sim- 
ulate. This can be done using sequences of collective 
pulses. Although the initial natural Hamiltonian restricts 
what operators can be obtained (see Appendix), various 
control sequences have been proposed to realize a broad 



2 



set of operators [19-21] starting from either the secu- 
lar dipolar Hamiltonian or an Ising interaction. These 
multiple pulse sequences cannot however modulate the 
coupling strengths, which is instead our goal. Evolu- 
tion under a magnetic field gradient imposes a modu- 
lation of the coupling strength. Consider for example 
that we want to modify the isotropic XY Hamiltonian, 
H XY = EijhjiSfS] + SfS]). To first order, evolu- 
tion under the propagator U(t, r) — e^ tHzT e~ HxYt e lHzT , 
where H z = ^ WiS* is obtained by a gradient, is equiv- 
alent to the evolution under the effective Hamiltonian 

H' XY = E« MOSf Sf + S?Sy) cos(^t) 
+(SfS] - SfS*) sm(SuijT)], 

where dcoij = Uj — Wj. The modulation can then be re- 
peated to obtain a total propagator Uq — J\ h U(th,Th). 
Given a desired target Hamiltonian H d = djj (SfSJ+ 
S?S]) + d-jiSfS] - SfSJ), we thus obtain a set of equa- 
tions in the unknowns {th,Th} that allow Hamiltonian 
engineering. To simplify the search for the correct tim- 
ings, we can further apply a filter that cancels all un- 
wanted couplings and use the equations above to only 
find the remaining, non-zero coupling strength. The fil- 
ter is obtained by creating a dynamic grating: instead of 
the propagator Uq, we evolve under N cycles (while re- 
ducing the times th th/N) with a gradient modulation 
U = rife=i e~* * U^e * T , which imposes a weighting 
factor Qij to the couplings, analogous to a dynamical im- 
plementation of a Bragg grating [7], 

We now make these ideas more concrete by consider- 
ing a specific example, the engineering of an Hamiltonian 
allowing perfect QIT in mixed-state spin chains, thus en- 
abling room temperature quantum communication [21- 
26]. For lossless transport, the simplest engineered n-spin 
chain consists of only NN couplings that vary paraboli- 
cally along the chain, dj = dWjJri — j) [18, 27]. This 
ensures perfect transport in a time T = n/(2d). How- 
ever, manufacturing chains with such restricted coupling 
topologies is a hard engineering challenge due to fabrica- 
tion constraints and the intrinsic presence of long-range 
interactions [28, 29]. On the other hand, regular, lattice- 
based spin networks (see e.g. Fig. 1) are found ubiq- 
uitously in nature. Our method can then be used to 
dynamically engineer the optimal Hamiltonian starting 
from a complex spin network. 

Filtered engineering for QIT - The target Hamilto- 
nian for QIT in a n-spin chain is H d — Y^jZi dj {S^ S^ +1 ~ 
SjSj +1 ). We consider a dipolar ly coupled spin network 
with Hamiltonian 

H = H int +H z =J2 b ij(3S?S?-S i -S j ) + J2"iS?, (!) 

ij i 
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FIG. 2. The filtered engineering control sequence consists 
of alternating blocks of evolution under the gradient T-i z for 
times Tj and DQ evolution Hbq of duration tj/N. The fun- 
damental cycle is repeated N times. 

where the spatial tagging of the frequencies is achieved 
by applying an appropriate magnetic field gradient. The 
double-quantum Hamiltonian Huq = J2i<j bij(Sf — 
SfSj) can be obtained from the secular dipolar Hamilto- 
nian via a well-known multiple-pulse sequence [19], which 
at the same time cancels the term H z . Importantly, the 
sequence allows time-reversal of the evolution by a simple 
phase shift of the pulses. We can further impose evolu- 
tion under the field gradient only, T-L z , by using homonu- 
clear decoupling sequences, such as WAHUHA [30, 31] 
or magic echo trains [32]. Alternatively, homonuclear 
decoupling can be avoided by shifting the DQ sequence 
off-resonance during the free evolution period (see Ap- 
pendix) . 

A control sequence achieving filtered Hamiltonian en- 
gineering consists of alternating blocks of evolution un- 
der % z and double quantum excitation %dq {mix- 
ing period, see Fig. 2). We analyze the net ef- 
fect of this pulse sequence using average Hamilto- 
nian theory [31]. Let us consider for simplicity a se- 
quence with only two mixing and free evolution blocks. 
Then setting U z (t) = exp(— irHz) and {7dq(£) = 
exp(— itHoQ), the propagator corresponding to N cycles 

is, U N = [U z ( n )Uv Q (ft) U z (t 2 )U U q (I) ] W , which can 
be rewritten as, 

U N = [UMUuq (ft) C/t( n) ] (2) 
x [U,(n + t 2 )Udq (ft) C/](n + r 2 )] 
x • • • x [U Z ((N - l)r)U m (ft) U$((N - 1)t)] 

where r = n + r 2 . Now, [U z (n)U D Q (ft) t/J(n)] = 
exp (-i^H m (Tt)), where H m (ri) = 

E t<J b tJ (S+S+e^ s >i + SrSje-^) is the toggling 
frame Hamiltonian with 5y = u)i + Uj . Employing the 
Suzuki- Trotter approximation [14], the propagator Un is 
equivalent to evolution under the average Hamiltonian, 

H = J2 ^+5+ (tie inS « + t 2 e^ +T ^) Gij + h.c. 
i<3 

where G l3 = e <(^-i)^/a ^^g) is the dynamical 
Bragg grating. 

In general, for a sequence consisting of free pe- 
riods t z = {r%,--- ,tl} and mixing periods t m = 
{<!,••• , ix,}, the average Hamiltonian is H = 
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FIG. 3. Left panels show control sequences for filtered en- 
gineering in a (a) 5 spin chain and (b) a 6 spin chain. The 
explicit values of t m are in . In the right panels, the phase ac- 
quired by the S^S^ term of the toggling frame Hamiltonians 
H mj (arising after each of the L mixing periods) is depicted 
on a unit circle [35]. To retain the DQ form of the average 
Hamiltonian, the phases have reflection symmetry about 2tt 
(see also ). 
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FIG. 4. Engineering filter amplitude \FijQij | for a 5-spin chain 
after different cycle numbers N. A single cycle creates the 
weighting function Fij shown in blue, which is transformed 
to sharp peaks at the ideal engineered couplings (red circles) 
as the number of cycles increases. 



T,i<jStSj'F ij (T z ,t m )g ij (T) +h.c. where r = J2f=i T j 
and we define the weighting function [7, 33], 

(t z , t m ) = j± tk exp (i5ij J2 k h= i Thj (3) 

k 

The grating Qij forms a sharp filter [7, 34] with maxima 
at r5ij — 2rmr. We assume to apply a linear ID-magnetic 
field gradient along a selected chain of spins in the larger 
network, such that the frequency of spin j is luj = jui—uiQ, 
where ujq is the excitation frequency. Each spin pair term 
of the DQ Hamiltonian acquires a spatial phase tag under 
the gradient field: if lot = tt and 2loqt = 37r — 2m7r the 
NN couplings are preserved, while the NNN couplings 
lie at the minima of the grating, and are decoupled (see 
Fig. 4). Other non-NN, off-chain couplings (see Fig. 1), 
lie at the grating side-lobes, and have greatly reduced 
amplitudes for large N, falling off linearly with ./V (see 
Appendix). Note that a suitable choice of gradient and 
free evolution could impose similar phase tags to other 
mixing Hamiltonians - for example, through a quadratic 
gradient for the isotropic XY Hamiltonian. 

Following the filter, the weighting function can 
be constructed to yield the ideal coupling strengths re- 
quired for perfect transport. Given the toggling frame 
Hamiltonians, H m {T^), we have a set of 2n equations 
(for an n-spin chain), 

Y,h=i sin(uj(2i + l)T h )thh,i+i = 0, Vi 

Ea=i c °s(w(2« + l)T h )t h b it i + i oc di, Vi, 

with 2L unknowns for L time steps. The number of con- 
ditions (and thus of time steps) can be reduced by ex- 
ploiting symmetry properties. For example, by imposing 
a gradient symmetric with respect to the center of the 
chain, it could be possible to automatically satisfy most 
of the conditions in Eq. 4 and only L = |~n/2] time steps 



would be required. Unfortunately this solution is practi- 
cal only for some chain lengths (see Appendix); we thus 
focus on a suboptimal but simpler solution. Let us first 
determine r z and t m for an odd n-spin chain. To enforce 
the mirror symmetry of the couplings dj [36] and en- 
sure that the average Hamiltonian remains in DQ form, 
we imposes mirror symmetry to the times, tj — th-ji 
while the gradient evolution intervals are Tj/t — 3/n for 
j = (L + l)/2 and tj/t = 1/n otherwise (see Fig. 3). 
This choice yields a simple linear system of equations for 
L = n — 2 mixing periods t m , 

An intuitive phasor representation [35] of how the evolu- 
tion periods exploit the symmetries involved is presented 
in Fig. 3 and in the Appendix. For even spin chains, 
the solution can be derived analogously with L = n + 1, 
<L = J2j dj and Tj = for j = L/2 + 1, and 
t z Jt — 1/n otherwise. 

When applied to a n-spin linear chain, the tuning ac- 
tion of Fjij + x\Qjtj + i\ is remarkably rapid, achieving per- 
fect transport fidelity in just a few cycles (Fig. 5). In- 
creasing the number of cycles reduces the error in the 
Trotter expansion improving Fjij + x\ (as shown in Fig. 
5. a) as well as improving the filter selectivity Gj(j+i) (see 
Fig. 5.b). In general, the peak width of the Bragg grat- 
ing decreases with the number of cycles as 2ir/N [7], im- 
proving the grating selectivity linearly with N (see Ap- 
pendix). Indeed, if all couplings are included, about n 
cycles are required for almost perfect decoupling of the 
unwanted interactions (fidelity > 0.95 in Fig. 5). 

The highly selective grating also avoids the need for 
the chain to be isolated and for the surrounding net- 
work to have a regular structure. Any spin that lies in 
between two nearest-neighbors of the chain can be ef- 
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FIG. 5. Minimum transport infidelity, with filtered engineer- 
ing as a function of cycle number N for a n-spin dipolar chain 
with (a) NN couplings only and (b) all couplings included. 
Here we calculated the fidelity F = Tr {USfU f S*} /2 n . Al- 
most perfect fidelity is achieved even for long chains with just 
a few filter cycles. 
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FIG. 6. (a) Transport fidelity of polarization in the network 
of Fig. 1 as a function of mixing time for N = 30. Almost 
perfect fidelity is reached for t = 7r/26j( i _|_ 1 ). (b) Variation of 
maximum fidelity with disorder in the surrounding network. 
The spins of the surrounding network are displaced by 8 ■ r, 
where S is the disorder strength, r is a random number in 
[—x/2,x/2] with x being the separation between NN chain 
spins, and we consider 30 random manifestations of the net- 
work. The filter improves as the number of cycles is increased. 



fectively decoupled, as shown in Fig. 6, where as small 
as N = 20 cycles are sufficient for robust QIT. However, 
in the presence of disorder in the couplings of the chain 
spins one needs to compromise between broader grating 
peaks (via small N) and poorer decoupling of unwanted 
interactions. 

Approximation mechanism - The control sequence is 
designed to engineer only the first order average Hamil- 
tonian H in Eq. (3). Higher order terms contribute with 
errors scaling as 0(tktk+i/N 2 ), as arising from the Trot- 
ter expansion [13, 14]. Consider for example the propa- 
gator for a 5-spin chain, 

U N = [e-^Me-^^^e-^^f, (5) 

where t\ — tt/(5uj) and — 2tt/(5uj) (Fig. 3). This 
yields the desired average Hamiltonian H with an error 
0{t\/N 2 ) for the first product, and 0{2t 1 t L /N 2 ) for the 
second. Increasing N improves the approximation, at the 
expense of larger overhead times (Fig. 5). In addition, 
the chosen scheme achieves remarkably good fidelities 
even for small N, since by construction, tj <C th ~ T/N. 
In essence, the system evolves under the unmodulated 
DQ Hamiltonian during tjj, yielding the average coupling 
strength, while the tj periods apply small corrections re- 
quired to reach the ideal couplings. 

Symmetrizing the recoupling sequence would lead to a 
more accurate average Hamiltonian because of vanishing 
higher orders [14, 37]. However, this comes at the cost 
of larger free evolution times t z and a larger number of 
periods of the unsymmetrized sequence could be fit into 
this overhead time. 

Experimental viability - We now consider the feasi- 
bility of experimentally implementing the filtered engi- 
neering sequence and show that high fidelity quantum 
transport at room temperature would be possible with 
current technology. 



For simplicity we assume that the spin chain or net- 
work is constructed in a physical lattice of NN separa- 

tion rn, leading to a NN coupling strength b = ^p--V- 



If an ideal n-spin chain can be fabricated in this lattice 
with maximum coupling strength b, the transport time 
required would be = ^ [22]. Alternatively, per- 
fect state transfer could be ensured in the weak-coupling 
regime [21, 38-40], with a transport time T wca k = 
where r>l ensures that the end-spins are weakly cou- 
pled to the bulk-spins. We compare and T wea k to the 
time required for N cycles of the engineering sequence, 
T ong . Since ij, 3> tj, to a good approximation the total 
mixing time is Ntj, < \Zj( n — j) l n ~ Since 
t = tt/lu, we have 



The overhead time — depends on the available gradi- 
ent strength uj. Since we can take Y ~ n for the weak 
regime [21] and N rj n for filtered engineering, a gra- 
dient larger than the NN dipolar coupling strength b 
would allow for faster transport in the latter case. As 
noted above, the chosen solution is not time-optimal and 
further speed-ups could be obtained for particular chain 
lengths. 

For concreteness, let us consider the crystal lattice of 
calcium Fluorapatite (FAp) [Ca 5 (P0 4 ) 3 F] [41, 42] that 
has been studied for quantum transport [43, 44]. It con- 
sists of 19 F nuclei that form linear chains along the c-axis, 
each surrounded by 6 other chains. The intra- nuclear 
spacing within a single chain is ro = 0.3442nm, which 
corresponds to b = 1.289kHz, while the inter-chain cou- 
pling is about 40 times weaker. Gradient strengths of 
comparable magnitude can be made by Maxwell field 
coils [45]: for instance, Ref. [46] demonstrated gradient 
strengths of 5.588 x 10 8 G/m over a 1mm 3 region; this cor- 
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responds to u = 0.7705kHz. Far stronger gradients are 
routinely used in magnetic resonance force microscopy 
(MRFM) [47]; for instance, dysprosium-based magnetic 
tips [48] yield gradients of 60G/nm, linear over distances 
exceeding 30nm, corresponding to cj = 82.73kHz. Con- 
servatively, one could estimate cj = 25kHz, for which 
usual NMR 7r/2-pulse widths of about 0.5/xs, would al- 
low sufficient bandwidth for coherent control on chains 
exceeding n = 50 spins. In order to ensure that the ex- 
perimental time does not exceed the decoherence time 
of the sample, one could apply homonuclear decoupling 
[31, 32] in the free evolution periods to reduce the effec- 
tive line width by at least an order of magnitude. For 
FAp, evolution under the DQ Hamiltonian survives for 
about 1.5ms [49] as limited by pulse errors; we anticipate 
that decoupling during the U z periods could increase this 
easily to 15ms [7, 32]. With u = 25kHz, and N = 30 cy- 
cles, nearly lossless transport should be possible for a 
25-spin chain. Note that in this case the free evolution 
period per cycle is 125/xs, allowing at least a few cycles 
of WAHUHA decoupling per U z evolution. 

Conclusion - We have described a method for quan- 
tum simulation without any local control, relying on the 
construction of time domain filter and weighting func- 
tions via evolution under a gradient field. The method 
was applied to engineer spin chains for perfect transport, 
isolating them from a large complicated networks. We 
showed that robust and high fidelity quantum transport 
can be driven in these engineered networks, with only 
experimental feasible control. 

Appendix 
Operator engineering 

Nuclear Magnetic Resonance has a long tradition of 
control sequences able to modify the naturally occurring 
Hamiltonians into desired operators; the most prominent 
application is in the refocusing of unwanted interactions, 
which has been further developed and is now a common 
technique in quantum information under the name of dy- 
namical decoupling. Here we extend these decoupling se- 
quences to more general two-body Hamiltonians as well 
as to the task of creating a desired operator. 

We consider general Hamiltonian for 2 spin-|: 

H = Si ■ <?i + w\ ■ a\ + a\ ■ d • <?2, (7) 

which can be rewritten in terms of spherical tensors Ti >m 
(see Table I): 

H = ^(-l) m A- m T *,m (8) 

1,171 

where the coefficients A;._ m depend on the type of spin- 
spin interaction and the external field. This notation 



is useful when considering rotations of the Hamiltonian. 
We consider only collective rotations -as easily achievable 
experimentally- that conserve the rank I. 

The goal of Hamiltonian engineering by multiple pulse 
sequence is to obtain a desired Hamiltonian from the nat- 
urally occurring one by piece- wise constant evolution un- 
der rotated versions of the natural Hamiltonian. We thus 
want to impose 

RkH-natRk = Udes, (9) 

k 

where Rk are collective rotations of all the spins, which 
achieves the desired operator to first order in a Trotter 
expansion. Rotations can be described by the Wigner 
matrices D l mn (R k ) as: 

R k Ti, m R\ = YSrVTEL,JRk)Ti, n (10) 
thus Eq. (9) can be written using the spherical tensors as 

k,n l,m l.m 

(ii) 

We thus obtain the set of equations 

YX-l) m Al? m DL,n(Rk) = (-l) n <-V (12) 

k,m 

We note that since collective pulses cannot change the 
rank I, there are limitations to which Hamiltonians can 
be engineered. In particular, T 0Q commutes with col- 
lective rotations: its contribution is thus a constant of 
the motion and, conversely, it cannot be introduced in 
the desired Hamiltonian if it is not present in the natural 
one. For example, an Ising Hamiltonian Hj = a z o~ z is ex- 
panded as Hi = (Too +v222o)/v3) so that only the sec- 
ond part can be modulated. Conversely, the secular dipo- 
lar Hamiltonian is given by T20, thus it cannot produce an 
Hamiltonian containing Too (although it can be used the 
create the DQ-Hamiltonian, Ruq = T22 -\-T2,—i). Group 
theory methods can be used to help solving Eq. (12) by 
reducing the number of conditions using symmetries. For 
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TABLE I. Spherical tensors for two spin-1/2 (a and b). a a 
are the usual Pauli operators. 



G 




FIG. 7. Phasor representation of the recoupling sequence for n — 5. The phase acquired by the toggling frame Hamiltonians 
H rnj after each mixing period in the sequence (a) is depicted on a unit circle for each of the NN terms Sj^u+i) in (b). The 
unit circle is traversed anti-clockwise from H mi to H mL . Mirror symmetry of the engineered couplings about the center of the 
chain is ensured since {Sf Sj, S^} and {S^ , S^} have identical phasor representations. The DQ form of the average 
Hamiltonian is ensured by forcing that the phasors for each of the Sj'Sf-^ terms have reflection symmetry about 2ir. 



Chain Length 


Mixing periods t m 


4 


{0,0.067,0.067,0,1.866} 


5 


{-0.201,-0.201,2.1242} 


6 


{0.0987, 0.1559, 0.0987, 0.0987, 0.1559, -0.0987, 2.6882} 


7 


{-0.369, -0.0652, -0.0652, -0.369, 2.8806} 


8 


{-0.2169, 0.1340, 0.2169, 0.1182, 0.1182, 0.2169, 0.1340, -0.2169, 3.4957} 


9 


{-0.5269, -0.1298, -0.0316, -0.0316, -0.1298, -0.5269, 3.6492} 


10 


{-0.3433, 0.0995, 0.2116, 0.2524, 0.1316, 0.1316, 0.2524, 0.2116, 0.0995, -0.3433, 4.2963} 


11 


{-0.6801, -0.1920, -0.0666, -0.0183, -0.0183, -0.0666, -0.1920, -0.6801, 4.4231} 



TABLE II. Mixing times t m in units of 1/d (as used in Fig. 5 and Fig. 6). Note that tL 2> tj for j 7^ L. 



example, the DQ Hamiltonian can be prepared from the 
secular dipolar Hamiltonian by using a simple sequence 
consisting of two time intervals, ii = £2/2 with the 



Hamiltonian rotated by R2 = f in second time period, 



to yield: T 2 nt x + 




2) — 5^0,0 



^2 CCHdQ- 



Symmetrized versions of this simple sequence are rou- 
tinely used in NMR experiments [50]. 



Phasor representation of the engineering sequence 

The filtered engineering sequence has an intuitive ge- 
ometric visualization in terms of phasors, which demon- 
strates the symmetries involved in the QIT Hamiltonian 
and the subsequent choice of free evolution periods r z . 
Let us consider for example the 5-spin sequence (Fig. 
7a), and the toggling frame Hamiltonians H m . during 
each mixing period. The phases acquired by the Hamil- 
tonians are represented on a unit circle for each of the 
Sj'Sj'^ terms (Fig. 7b). The phasor representation 
makes clear the symmetries involved in the filter en- 
gineering. Since the ideal couplings dj are mirror sym- 
metric about the center of the chain, the construction of 
the filter i^Yj+i) should match this symmetry. This can 
be ensured in two ways. 

We could select a gradient centered around the mid- 



dle spin in the chain. Then uij = —uj n -j, and the set of 
equations in (4) reduces to only n — 1 equations (instead 
of 2(n - 1)) if t k = t L - k and n m (r k ) = H m (T L _ fe ) T . 
This second condition further ensures that the first set of 
equations in Eq. (4) are satisfied (the condition is easily 
satisfied by setting t z k = —t z L _ k with t z k = r k — T^-i and 
t\ = Ti). Thus we only have \(n — l)/2] equations to 
be satisfied and correspondingly only \(n — l)/2] mix- 
ing periods. Although solutions can be found with this 
scheme, only for n — 4, 5 all the times are real. For 
larger spin chains it is possible to find all real solutions 
for L = \(n — l)/2] + 2 mixing periods, although the 
system of equations quickly becomes intractable. 



A second strategy is to use a gradient with minimum 
at the first spin in the chain, L = n — 2 mixing periods 
(for n odd) and choose Tj/r = 3/n for j = (L+ l)/2 and 
Tj/t = l/n, and tj = t( L _jy Moreover, since we want to 
retain the DQ form of the average Hamiltonian H, the 
toggling frame Hamiltonians H m (r k ) and H m (TL_ k ) are 
such that their phasors have reflection symmetry about 
2-7T for each Sj~Sj~ +1 term. This ensures that only the 
(Sj~S^ +1 + SJ SJ +1 ) cos[t<5j( :(+1 )] term survives, leading 
to an effective DQ form. 
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FIG. 8. (a) The amplitude of the dynamic Bragg grating, 
\Gij\ for N = 5 cycles. There are iV — 1 minima at 2jn/N, 
and N — 2 secondary maxima (sidelobes). The width of the 
primary maximum is 2n/N, and decreases with the number 
of cycles making the grating more selective, (b) Right panel 
characterizes the decoupling efficiency in terms of the mean 
amplitude of grating \Qij\ in the primary maxima (i.e the 
range [0,7r/iV]U[2(JV-l)7r/iV, 2-k/N]), and the side lobes (the 
range [ir/N,2(N — 1)tt/N]). Increasing the number of cycles 
decreases the relative sideband power, and hence increases the 
grating selectivity. 



Characterizing grating selectivity 
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FIG. 9. A possible implementation of the U z blocks is 
by shifting the offset frequency of the DQ excitation cjn to 
■mrmbjtj+iy Secular truncation ensures that even an offset by 
a few filter periods is sufficient to obtain fidelities (solid lines) 
comparable to when the U z blocks are created by dipolar de- 
coupling (dashed lines). We considered in (a) the network in 
Fig. I and in (b) a 9-spin chain with all couplings. 



vantage of enhancing the net decoherence time of the 
system, allowing for transport in longer engineered spin 
chains. 



Let us now characterize the selectivity of TV-cycles 
of the Bragg grating = J2k=o e ikT5i i . Fig. 8(a) 
shows the grating amplitude \Gij\ for N = 5 cycles. 
In general, the width of the primary maximum, 2-k/N, 
decreases with N and this increases the selectivity of 
the grating. Concurrently, the relative amplitude of 
the N — 2 sidelobes decreases with N. Specifically, 
consider Fig. 8(b), which compares the mean grating 
amplitude in the primary maxima i.e. in the range 
[0,tt/7V] u [2{N - 1)tt/N,2tt/N], and the sidelobes in 
the range [w/N,2(N — 1)tt/N]. From the properties of 
the Bragg grating, we find that the decoupling efficiency, 
characterized by the relative amplitude of the amplitude 
in the primary maxima with respect to the sidelobes, in- 
creases linearly with N. 



Off resonance DQ excitation to construct U z 

In the main text, we assumed that the free evolu- 
tion propagators U z were constructed by employing a 
homonuclear decoupling sequence (for eg. WAHUHA) 
during the t z intervals to refocus the internal Hamilto- 
nian couplings, while leaving the action of the gradient. 
An alternate method is to use the DQ excitation sequence 
even during these periods, but shifting the offset fre- 
quency wo far off-resonance. This is possible since the 
recoupling filter is periodic with a period 2mir thus we 
can shift ujq while retaining identical filter characteris- 
tics. Fig. 9 shows that an offset a few times larger than 
the coupling is sufficient to reach almost free evolution. 
Still, employing a decoupling sequence requires almost 
the same control requirements, and has the added ad- 
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